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Abstract 

The electromagnetic response of a nanoparticle of an ion-doped polymeric elastic 
insulator, commonly called as an electret, is considered in the continuum model of a 
uniformly charged elastic sphere. The spectral formulae for the frequency of optically 
induced spheroidal and torsional shear oscillations driven by bulk force of elastic and 
dielectric stresses are obtained in analytic form. Particular attention is given to re- 
laxation dielectric mode of the electrostriction response and its stability in the lowest 
quadrupole mode. The practical usefulness of ultrafine particles of electrets as a bio- 
labels capable of accumulating likely-charged inclusions uniformly dispersed over the 
spherical volume of an elastic matrix is briefly discussed. 
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1 Introduction 

The phenomenological continuum-mechanical models of vibrating liquid drop and solid globe 
provide proper account of gross features of resonance-like excitations in electromagnetic re- 
sponse of ultra fine particles. The spectra of these excitations are currently measured by 
methods of Raman scattering and inelastically scattered neutrons. The physical idea un- 
derlying these models is to identify the frequency of ac field incident on a nanoparticle with 
frequency of material oscillations in the particle regarded as a spherical piece of either liq- 
uid or solid continuous medium. A canonical example is fluid-mechanically computed Mie's 
spectral formula for the frequency of surface plasmons in nanoparticles of highly conducting 
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metals like silver and gold [1-5]. Another example is the solid-mechanically calculated fre- 
quency spectra of optically induced elastic oscillations, acoustic phonons, in a nanoparticle 
of a dielectric solid [6-8]. 

In this work we consider a model of resonant response of a nanoparticle of charged elec- 
tret (insulating elastic matrix doped by likely-charged ions homogeneously distributed over 
spherical volume) responding to electromagnetic load by shear electro-mechanical fluctua- 
tions which are described in the approximation of a uniformly charged elastic continuum. 
To illuminate the difference between the electric reactive forces driven by optically-induced 
electro-mechanical vibrations in a nanoparticle of a uniformly charged electret (dielectric 
material) and plasma oscillations of free electrons against immobilized ions in a nanoparticle 
of noble metals (highly conducting materials), in Sec. 2 a brief outline is given of classical 
electron theory of the surface plasmons and solid-mechanical theory of acoustic phonons. In 
Sec. 3, the frequency spectra of electrostriction spheroidal and torsional vibrational modes 
in the ion-doped dielectric nanoparticle of electret are obtained followed by discussion of 
electro-elastic stability of such a nanoparticle to optically induced distortions, the issue of 
crucial importance to their use as biolabels. The results are briefly summarized in Sec.4. 

2 Resonance frequencies of surface plasmons and acous- 
tic phonons 

Within the framework of the classical electron theory of metals, the low-frequency elec- 
tromagnetic response of metallic nanoparticle to long-wavelength ac field is modeled by 
oscillations of electron liquid against crystal lattice treated as a static uniformly charged 
ionic background and described by coupled equations of electron fluid-mechanics and ionic 
electrostatics 

STT-. * /-,\ 

p-^^pedE, p^m^Ue, Pe^-eUe- (1) 

The emergence of Coulomb restoring force is attributed to a small amplitude fluctuations 
of intrinsic electric fleld 5E (between free electrons and ions) superimposed on electrostatic 
fleld of ionic jellium, Eq which is computed from the equation V • Eq = 47rpi, where pi — enf, 
the particle number of electrons Ug equals to that of ions rii due to electro-neutrality of metal 
nanoparticle. The effective mass of electron m* is unique to each speciflc metal. The effect 
of resonance photo- absorption emerges when the frequency of incident on particle ac fleld 
SFi becomes equal to the frequency of oscillations of the velocity Sv of electron flow obeying 
the vector Laplace equation: 

V^5E = 0, V^5v = 0. (2) 
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JP^dV = J Pe(5v-5E)rfV. (3) 



These equations can be thought of as long wavelength limit of the Helmholtz equation for 
standing waves and describing electromagnetic response by long wavelength vibrations of 
electric field and electron flow, Rayleigh regime. In this regime the frequency of oscillations 
of free electrons against ionic background can be computed by the energy method which is 
based on the equation of energy balance 

d f p6\-'^ 
di 

The frequency of surface plasmons as a function of multipole degree £ is given by the well- 
known Mie spectral formula [1-4]: 

47rnee^ 
'2£+l ~P ml 

The dipole surface plasmon resonance whose energy is given by, Euiei^^ = 1) = fv^p/V^-, is 
well singled out in the cross sections of photo-absorption by silver and gold nanoparticles. 

Also relevant to our further discussion is the continuum, solid- mechanical, description 
of acoustic phonons, quanta of low-frequency elastic oscillations of solid. The mathematical 
treatment of such vibrations rests on canonical equation of solid-mechanics for the field of 
material distortions Ui{r,t) 

PUi = "^kCTik, CTik = '^■pUik, Uik = ^[VjMfc + VfcMi], Ukk = ^kUk = 0. (5) 

where p stands for the density the shear modulus p relates applied shear stresses aik and 
resulting shear or strains Uik in accord with the law of Hooke. The conservation of energy is 
controlled by equation 

d f pii^ f 1 

^ / (TikiiikdV, Uik = -[Viiik + VkUi]- (6) 

The fluctuating fields of material displacements Mj(r, t) and strains Uik{r,t) can be conve- 
niently represented in the following separable form 

u(r, t) = a(r) a(t), Uik{r, t) = aik{r)a{t), a^k = ^[Vi + a^] (7) 

with help of which the equation of energy balance (E]) is reduced to equation for a{t) having 
the form of standard equation of linear oscillations: 

Ma(t) + = 0, M = j pai{r)ai{r)dV, K = 2 j paikir) aik{r)dV. (8) 

In Rayleigh's regime, the vibrational response of an elastic sphere is characterized by two 
fundamental modes of shear elastic oscillations of the nodeless fields of material displacements 
which are specified by two fundamental solutions of the vector Laplace equation 

V2u(r,t) = ^ VVr) = 0. (9) 
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The instantaneous displacements in spheroidal mode of irrotational oscillations - is given 
by poloidal field and in torsional mode of differentially rotational oscillations - by the 
toroidal fields, respectively [7j: 

a, = VxVx(rx), a^ = V x (rx), xW = AT^ / ^^(cos 0). (10) 

Henceforth Pe{cos6) stands for the Legendre polynomial of multipole degree £ and A/^ is the 
arbitrary constant. The spectral equations for the frequencies of the even parity spheroidal 
Us and the odd parity torsional ujt modes of acoustic phonons are given by [6l [7] : 

Us{e) = a;o[2(2£ + !){£- 1)]'/^ = Uoim + 3)(£ - 1)]'/^ i > 2, (11) 

a;o = |, Q = y^ (12) 

where q is the speed of elastic shear wave in the material bulk and R stands for the radius 
of nanoparticle. These last spectral equations are central to spectroscopic analysis of optical 
modes excited by inelastically scattered x-rays and neutrons. The purpose of the above 
outlined models is to illuminate similarity and difference between spectral equations for the 
frequencies of surface plasmons, equation ([5]), and the frequency spectra of electrostriction 
modes in a nanoparticle of a uniformly charged electret which is the main subject of the 
next section. 



3 Electrostriction mode in optical response of nanopar- 
ticle of a uniformly charged electret 

The optically-induced, by ac electromagnetic field, electro-mechanical distortions in the vol- 
ume of uniformly charged electret are characterized by the tensor of dielectric stresses (e.g., 

my- 

Pik = ^ [Ei 5Dk + Ek 5D, - {E, 5Dj)5,k] (13) 

where Ei are components of electrostatic field produced in the particle volume by extraneous 
charge uniformly distributed with the charge density pe- 

V-E = 47r— , VxE = (14) 
e 

whose solution is well known: 

E(r) = -V$(r), $(r) = -|^pe(r2-3i?2), [E, = l![^r, i?, = 0, i?^ = 0]. (15) 

o6 o 6 
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The electro-mechanical effect is described by constituting equatioiu 

6Di = -inpeUi (16) 

showing that optically induced ffuctuations of dielectric induction 6Di are linearly propor- 
tional to material displacements Ui and inextricably related to the storage of extraneous 
charge uniformly dispersed with density pe over the sample volume; note, the dielectric 
materials can accommodate only extraneous charge. 



3.1 Dielectric relaxation mode 

We confine our analysis to the Rayleigh's regime of optical perturbation resulting in non- 
compressional fluctuations of the electret material (the charge density remains unchanged 
6pe = —pe Vfc Uk = 0) which are described by nodeless field of material displacements Ui 
obeying the vector Laplace equation, V^u = 0. The purpose of this subsection is to elucidate 
the difference between plasma oscillations of free electrons against ionic jellium in a metal 
nanoparticle which are restored by electrostatic force of the form fi = peEi, Eq.([l]), and 
the shear oscillations in a nanoparticle of uniformly charged electret which are restored by 
reactive force of dielectric stresses // = Vfc Pik which are governed by equations 

d f pu^ f 1 

pui = VkPik, ^ / — dV = - Pij iiij ciV, iiik = - [Vj Uk + Vfc iii] , Ukk = 0. (17) 

The rightmost equation shows that optically induced dielectric stresses pik are accompanied 
by shear deformations Uik- The eigenfrequency can be computed with help of the above 
expounded Rayleigh's energy method whose key point is the separable representation of 
fluctuating material displacements and strains in the form given by equations ([7]). Then, for 
the perturbation-induced dielectric induction one has 

5Di{r,t) = -A7rpeai{r)a{t) (18) 

and the tensor of dielectric stress is given by 



Pik{r,t) = [Tik{r) - ^Tjj6ik]a{t), Tik{r) = -pe[^i(r) afc(r) + Ek{r) ai{r)]. (19) 



^ The constitutive equation for SDi is compatible with the Maxwell equation V x = (47r/c)(5j + (l/c)(5D. 
Applying to this latter equation operator of divergence and taking into account that 6j = PeSv and the 
continuity equation of the charge conservation Spe = — V^j, we obtain SD — —iirpeSv. Bearing in mind 
that 6v = u and eliminating in the last equation the time derivative one arrives at ([16]). 
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Substitution ([7]), (fTSj) and (fT9l) into equation of energy balance in ( !T7|) leads to 

_ = ff=_ + ^ ^ ci(i)+^VW = 0, (20) 

M = fpa,a4V, /<-. = /r..(r)a..(r)dV (21) 

x;,fc(r) = -pe[-Ei(r) aA:(r) + Efc(r) ai(r)], a^^ = - [Vj + Vjfc aj] . 

The fact that such response is accompanied by internal shear deformations suggests that 
dielectric modes in question can be specified in a manner of vibrational modes in an elastic 
sphere, that is, as spheroidal and torsional ones. 

Dielectric spheroidal mode. In this positive parity mode the displacements are described 
by the poloidal (polar) vector field 

a, = TV; V X V X r / (cos 0) . (22) 
For the inertia, stiffness and frequency ujds of spheroidal dielectric mode we obtain 
/lf^-4vr.AA2ff2m^(i+l)! j,s_ 327rVg ,,^, £(£+l)(£^-l) 

^iW = -1-^(^-1), -^ = -^ (23) 
6 e p 

where Ud is the natural unit of dielectric fluctuations. 

Dielectric toroidal mode. The material displacements in negative parity torsional mode 
are described by the axial toroidal (axial) vector field 

at = MiV X r/P£(cos0). (24) 

Computation of the the inertia, stiffness and frequency Udt of torsional dielectric mode yields 

M' = A^pUfR^'^ \ ^^^ + ^) K\ = -^-^^UfR^'-'^ ^(^'-1) 



(2£ + l)(2£ + 3)' " ?> e ' (2£ + l)(2^ + 3)' 

^^.W = # = -^-J(^-l), ^^ = -^- (25) 
Mt 6 e p 

The obtained spectral equations represent one of the main newly obtained result of the 
presented theory showing that optical response of nanoparticle of electret are characterized 
by two different in parity modes, even-parity spheroidal dielectric mode and odd-parity 
torsion one. The basic dielectric frequency Ud depends upon the dielectric constant e in such 
a way that in the limit e — > oo, as is the case of conductors, Ud 0. This means that. 
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contrary to the surface plasmons in a nanoparticle of a highly conducting noble metals, the 
electrostatic fluctuations in a nanoparticle of a uniformly charged electret are manifested as 
relaxation modes, not oscillatory. It is relevant to note that the somewhat different method 
of computing the frequency of nodeless electro-elastic vibrations of a uniformly charged 
dielectric globe that has been developed in [101 [H] leads to identical analytic estimates for 
the frequencies of considered dielectric modes. The negative sign of squared frequencies, 
highlighting relaxation nature of dielectric mode, indicates that nanoparticle is unstable to 
optically induced electro- mechanical distortions^. In the next section we consider oscillatory 
response of such particle which is governed by superposition of two restoring forces, namely, 
the constructive elastic restoring force whose strength is dominated by gradient in shear 
modulus 

= VfeCTifc = 2uikVkfi{r) + 2fiVk Uik ^ 2MifeVfe/i, Uik = V'^Ui = 0. (26) 

and the above destructive force of dielectric stresses, // = VkPik and and derive condition 
of stability which imposes severe constrain on the total charge and radius of the uniformly 
charged nanoparticle of electret which is of crucial importance to their practically useful 
applications. 

3.2 Frequencies of electrostriction modes 

The optically induced vibrations driven by combined the above forces of gradient-shear 
elastic stress /" and dielectric stresses /' are governed the following equations of electro- 



■^The origin of this instability can be traced from equations of electro-elastic solid mechanics describing 
propagation of perturbations in bulk of uniformly charged electret with constant shear modulus /i which is 
given by 

pill ^ ^kCTik + Vi pife = M u + (47r/e) pi u. 

For the plane- wave form of externally induced displacements, u = expz(kr — ojt), the last equation leads 
to dispersion equation 



„2 



2 47r/9^ 

p e p 



where Ct is the speed of transverse shear wave in the bulk electret and ujd stands for the natural unit of 
frequency of dielectric oscillations. It follows that there is critical wave number (and corresponding 
wavelength Ac = 27r/fcc such that the electret material becomes electro-mechanically unstable when 



Understandably that this instability owe its origin to the presence of extraneous charges. 
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elastic dynamics and energy balance 



piii^ WkCTik + ^iPik = ^UikS/kfi + —Vk[Ei 5Dk + Ek (^A = -47rpe Wj, (27) 

^ j ^dV = -2 y iVkPi)uikUi(n^ - J Pe[Ei Uk + Ek Ui\ iiikdV (28) 
and making use of substitution tii(r, t) = aiiv) a{t) we obtain 

M. + i.«^0. M^lpa.a.^,^ K = K, + K.. (29) 

Ki^^2 J (Vfc//) Oife dV, Kd^-^ J Pe[Ei ttk + -Efc a,] [Vj Ofe + Vfe a^] cH^. (30) 
The mass parameters in spheroidal and torsional modes are given by 



and for the stiffness one has 

R 





R 



3 e 
167r2 p2 



(2£+l) 

?2-l) 



e \f2 o2e+3 ' 

3 e * (2£+l)(2£ + 3)" 



(31) 



(32) 
(33) 



These equations provide a basis for computing frequency 



M 



of electrostriction vibrational response of a nanoparticles of uniformly charged electret with 
the non-homogcncous shear modulus profile. 

As a representative example let us consider a case of nanoparticle with shear modulus 
profile given by 



(34) 



In this case the frequency of spheroidal electrostriction mode can be conveniently represented 
in the following form 



00. 



ul2{2£+\){£-\)r^ 
8 



3 
1 - 



3(2£ + 1) 



/5 



(35) 



and for the torsional electrostriction mode as follows 



,(2£ + 3)(£- 1) 

2(£+l) 
2 (2 £ + 3)(£- 1 
2(£+l) 



3 
1 - 



- 1) 



2(£+l) 
3(2£ + 3) 



/5 



(36) 



One sees that the lowest overtone of both spheroidal and torsional modes is of quadrupole 
degree ^ = 2. The absence of monopole £ = 0, breathing overtone, is the consequence of 
adapted approximation of incompressible matter. The dipole fields of both poloidal and 
toroidal displacements describe center-of-mass translation and rigid-body rotation, respec- 
tively, that is, the non-vibrational reaction of nanoparticle. 




Figure 1: Ratio of the squared frequency of spheroidal (upper panel) and torsional (down panel) 
electrostriction modes to the squared frequency of elastic shear oscillations as a function of multipole 
degree at indicated values of parameter of stability (3. 

In Fig.l we plot the ratio ojI/ujI and ojI/ojI as functions of the multipole degree of 
spheroidal and torsional vibrations respectively, showing that the larger the C. the higher is 
the frequency. Also, this figure exhibits strong dependence of frequencies upon the parameter 



Q = peV, V = (47r/3)i?l 



(37) 
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carrying information about total charge Q accumulated by particle of radius R and shows 
that the larger (5 (the large the ratio Q/ R^) the lower is the frequency. The most conspicuous 
feature of electrostriction response in question is that the lowest, quadrupole, overtones 
become unstable when the parameter (3 attains critical value [3 = Pc- Specifically, the 
spheroidal electrostriction vibrational mode becomes unstable, meaning uJs{^ = 2) = 0, 
when 



1 - 



3(2£+ 1 



-P 







£=2 



15 

Y' 



(38) 



The lowest quadrupole torsional electrostriction vibrational mode unstable, u!s{i = 2) = 0, 
when 



2(^+1) • 
3{2i + 3f_ 







Pi 



1=2 



7 
2' 



(39) 
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Figure 2: The frequency of quadrupole spheroidal and torsional electrostriction modes as a function 
of the stability parameter (3. The circles on OX-axis mark values of (3 at which the electrostriction 
mode becomes unstable undergoing transition from regime of oscillations to the relaxation regime. 

In Fig. 2 we plot the ratio of squared frequencies of quadrupole, i = 2, overtones of both 
spheroidal and torsion electrostriction modes as functions of the stability parameter /3 by 
highlighting the above critical values of stability parameter (3^ by circles on the OX-axis. In 
this points the electrostriction modes undergo transition from the oscillatory regime to the 
relaxation regime. This leads us to conclude that nanoparticle of uniformly charged electret 
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is stable to optically induced deformation oscillations when 

< C(v^) (40) 

where constant C falls in the range 5 < C < 10; understandably that the lowest /9c should 
be regarded as genuine critical value of this parameter. The practical usefulness of the above 
established conditions of instability is that it imposes severe constrain on the size and total 
charge accommodated by the electret nanoparticles and this must be taken into account in 
the process for technological fabrication. As a representative example, for the nano and micro 
dimensions 10~^ < R < 10~^ m whose total charge is of the order of Q ~ lO^e ~ 10~^^ (in SI 
units) from polymers with dielectric constant 1 < e < 15, and shear modulus 10^ < fi < 10^ 
(in SI units), from above obtaned condition it follows that 10~^ < Q/R^ < 1. So that for 
ultra fine particles of nano sizes this condition is always fulfilled. 



4 Summary 

In this, somewhat technical in mathematical details, paper the theory of electrostriction 
response of nanoparticlc of non-conducting clastic polymers capable of accumulating an ex- 
traneous charge to ac electromagnetic load has been developed. The prime purpose was to 
elucidate stability of such nanoparticles to optically induced distortions. This issue is of cru- 
cial importance to utilization of uniformly charged polymeric ultra fine particles (which are 
currently produced by jammed technique) as biolabels as well as for other hi-tech application. 
Based on general equations of electro-solid-mcchanics appropriate for an isotropic electrets 
capable of accumulating likely- charged inclusions uniformly dispersed over the volume, the 
model of spherical particle of such a material undergoing optically induced spheroidal and 
torsional oscillations restored by superposition of two forces has been considered, namely, 
the constructive elastic force dominated by gradient of shear modulus and destructive force 
of dielectric stresses has been considered. It is that, contrary to nanoparticle to noble metals 
(responding to ac field by surface plasmon resonances which are independent upon particle 
size), the optical response of nanoparticles of uniformly charged electrets strongly depends 
upon their size and total accumulated charge. The conditions of electro-elastic stability of op- 
tically induced electrostriction spheroidal and torsional vibrational modes has been obtained 
showing that both modes of such response undergo transition from regime of oscillations to 
the relaxation regime. This condition highlights the most important difference between op- 
tically induced electrostatic fiuctuations in the dielectric nanoparticle of electret which are 
of volume character and electrostatic, plasma, fiuctuations in the metal nanoparticle which 
are of surface character: in this latter case these oscillations are stable whereas in the former 
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case the electrostatic fluctuation can lead to disintegration of the electret nanoparticle by 
means of electric discharge. All these suggest if the condition of electro-elastic stability is 
fulfilled the optically induced electrostriction modes in question can be detected as resonant 
modes of photo absorbtion and the lowest overtones of these modes are of quadrupole degree. 
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